Validity of single-channel model for a spin-orbit coupled atomic Fermi gas near 

Feshbach resonances 
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We theoretically investigate a Rashba spin-orbit coupled Fermi gas near Feshbach resonances, 
by using mean-field theory and a two-channel model that takes into account explicitly Feshbach 
molecules in the close channel. In the absence of spin-orbit coupling, when the channel coupling g 
between the closed and open channels is strong, it is widely accepted that the two-channel model 
is equivalent to a single-channel model that excludes Feshbach molecules. This is the so-called 
broad resonance limit, which is well-satisfied by ultracold atomic Fermi gases of ^Li atoms and ''"K 
atoms in current experiments. Here, with Rashba spin-orbit coupling we find that the condition for 
equivalence becomes much more stringent. As a result, the single-channel model may already be 
insufhcient to describe properly an atomic Fermi gas of ^''K atoms at a moderate spin-orbit coupling. 
We determine a characteristic channel coupling strength gc as a function of the spin-orbit coupling 
strength, above which the single-channel and two-channel models are approximately equivalent. We 
also find that for narrow resonance with small channel coupling, the pairing gap and molecular 
fraction is strongly suppressed by SO coupling. Our results can be readily tested in atoms by 
using optical molecular spectroscopy. 



China 



I. INTRODUCTION 

As a realization of non-abelian gauge fields in neutral 
cold atoms [H-Q, spin-orbit (SO) coupled atomic gases 
have attracted a lot of attentions in recent years. The 
SO coupled bosonic gas of *^Rb atoms was first achieved 
by Spielman's group at National Institute of Standards 
and Technology (NIST) in early 2011 py]. The SO cou- 
pled atomic Fermi gas has also been realized most re- 
cently at Shanxi University f3| and at Massachusetts In- 
stitute of Technology (MIT) [5] with ''°K and '^Li atoms, 
respectively. These novel atomic gases have many in- 
teresting properties inherent to spin-orbit coupling, and 
have potential applications in future quantum technol- 
ogy. A well-known example is the emulation of the lonj 
sought topological superfluids and Majorana fermions 

which lie at the heart of_topological quantum infor- 
mation and computation (oi. flO]. 

Most of previous theoretical studies on SO coupled 
Fermi gases are based on a single-channel model [6, 11,- 
IbH - In this model, the interaction between atoms is de- 
scribed by a single parameter, i.e., the s-wave scattering 
length as- The scattering length can experimentally be 
tuned by using Feshbach resonances. As a result, in the 
absence of SO coupling, the Fermi gas can cross from 
a Bose-Einstein condensate (BEG) over to a Bardeen- 
Gooper-Schrieffer (BCS) superfluid [32], when the scat- 
tering length changes from positive to negative values. 
With SO coupling, the picture of BEG-BGS crossover 
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may be qualitatively altered. For example, for a partic- 
ular Rashba-type SO coupling, a new two-body bound 
state - referred to as rashbon - is formed [ll| - [l4j . By 
increasing the SO coupling strength, the system may 
change from a BGS superfluid to a BEG of rashbons, 
even on the BGS side with a negative scattering as < 
[l^ - [l^ . The pairing gap of this system will be signifi- 
cantly enhanced due to the increased density of state at 
the Fermi surface [l3|. An anisotropic superfluid due to 
the Rashba SO coupling has also been predicted [13.] . 

A more realistic and complete description of ultracold 
atomic Fermi gases near Feshbach resonances, however, 
should be the two-channel model, which includes both 
atoms in the open channel and Feshbach molecules in 
the closed channel In this model, in addition to the 
background s-wave scattering length abg between atoms, 
two other parameters are used in order to fully describe 
the interaction. These are the detuning energy of Fes- 
hbach molecules v and the channel coupling strength g 
between molecules and Fermi atoms. Therefore, the in- 
teraction of the system consists of two parts. The non- 
resonant part is the contact interaction between atoms 
with the strength determined by the background scatter- 
ing length, while the resonant interaction is induced by 
the coupling between molecules and atoms. Near Fesh- 
bach resonances without SO coupling, it is known that 
the single-channel and two-channel models are essentially 
equivalent when the channel coupling strength g is large 
enough fsS, l35|. This is the so-called broad resonance 
condition, satisfied by the Fermi gases of "^^K and ^Li 
atoms, which are so far the two main systems used in the 
cold-atom laboratory. 

In this paper, we aim to examine the equivalence of 
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the single-channel and two-channel models for a Rashba 
SO coupled Fermi gas near Feshbach resonances. This 
is by no means obvious, as fermionic pairing is notably 
affected by SO coupling at the BEC-BCS crossover. We 
use mean-field theory and focus on the most interesting 
resonant limit. Our results show that in the presence 
of SO coupling, the broad resonance condition is much 
more difficult to achieve. As a result, for an ultracold 
atomic Fermi gas of ''^^K atoms, which is known to be 
well described by the single- channel model without SO 
coupling, we may have to use a two-channel model al- 
ready at a moderate SO coupling strength. 

Our paper is organized as follows. In the next section 
(Sec. II), we introduce the model Hamiltonian. In Sec. 

III, we diagonalize the Hamiltonian by using mean-field 
theory to obtain the grand thermodynamic potential and 
solve the resulting coupled mean-field equations. In Sec. 

IV, we discuss the equivalence between the single-channel 
and two-channel models. In Sec. V, we show how to test 
experimentally the difference between the two models, 
by using optical molecular spectroscopy [36| . Finally, we 
summarize in Sec. VI. 



II. MODEL HAMILTONIAN 

We consider a three-dimensional (3D) resonantly- 
interacting atomic Fermi gas with Rashba-type SO cou- 
pling, described by the two-channel model Hamiltonian, 



where aj^^ is the creation operator for atoms with mo- 
mentum k in the spin state cr, Ck = fi.^k^/(2m) and 
</Jk = a'rg(A;2; + iky). To diagonalize this single-particle 
Hamiltonian, we introduce the field operators in the he- 
licity basis labeled by "±", 
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(4) 



with which the single-particle Rashba Hamiltonian be- 
comes diagonal. 



Hso = 



(5) 



Note that, in the helicity basis, the single-particle dis- 
persion relation now breaks into two branches: ek+ ~ 
Ek + h'^Xk±/{2m) for the upper branch and Ck- = Ck — 
h'^Xk±/{2m) for the lower branch. To describe Feshbach 
molecules, we use annihilation operators bq. The energy 
of molecules is denoted as 2]^, which after renormaliza- 
tion [ss'l is related to the detuning energy from threshold 
of Feshbach resonance Bq, i.e., 2i^o = Afj,{B — Bq), where 
Afi = 2fia — fJ-m is the magnetic moment difference be- 
tween the atomic {2fia) and bound molecular state (fim) 
(sH - The Hamiltonian of Feshbach molecules may be 
written as, 



H=HsO+Hrn+Hl, 



(1) 



where Hso and "Hm stand for the non-interacting Hamil- 
tonian of SO coupled atoms in the open channel and 
of Feshbach molecules in the closed channel, respec- 
tively. The interaction Hamiltonian Hi = T-Lam + Haa 
includes both the atom-molecule coupling between the 
two channels {Ham) and the background interaction be- 
tween open-channel atoms {Haa)- 

For atoms, we take the following single-particle Rashba 
SO Hamiltonian, 



n,n = 2iyY,blbq. 



(6) 



Finally, the interaction Hamiltonian is given by 'Hi = 

Haa +Ham, whcrC 



Haa - Ubg ^ a|j/2+kt"q/2-ki«q/2-k'iaq/2+k't (7) 
kk'q 
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Hso = 



2m 



(2) 



where kj^ = {kx,ky) and a±^ = {(7x,<^y) are respectively 
the in-plane momentum and in-plane Pauli matrix, and 
A is the Rashba SO coupling strength. Note that, the 
standard representation of the Rashba SO coupling is 
given by \{kyax — k^dy) [l3| . Here, for convenience we 
have performed a spin-rotation to rewrite the Rashba 
term into a slightly different but fully equivalent form 
X{k ax + kyUy) In the second quantized form. 



Hso = y^£kQ^ 

kcr 



kcr 



flko 



2m 



H.c. 



(3) 



is the non-resonant interaction between atoms, with 
strength given by the background s-wave scattering 
length after renormalization [33], Ubg — ATrh^abg/m, and 

Ham = 5 X! K«q/2+ktaq/2-ki + H.c] (8) 

kq 

is the resonant interaction between atoms and molecules, 
with strength parameterized by g. After renormaliza- 
tion, the magnitude of the channel coupling strength g is 
related to the width of the Feshbach resonance W, i.e., 
g = y/ AfiWUbg- We note that, in the two-channel model 
one may define an effective s-wave length [35], 



W \ g" m 
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III. MEAN FIELD THEORY 

We use the standard mean-field theory to solve the 
two-channel model Eq. ([T]), by assuming that all the 
molecules and Cooper pairs condense into the zero- 
momentum state. Thus, we set q = in the interaction 
Hamiltonian "H/. Here, we have excluded the possibil- 
ity of an inhomogeneous superfluid phase (i.e., q 7^ 0), 
which may exist in the presence of an in-plane Zeeman- 
field l'37'l. This is consistent with the two-body calcula- 
tion [13, [ill that the ground state of two particles in 
our Hamiltonian always has zero center-of-mass momen- 
tum. Following the procedure in Ref. (33| . we introduce 
the following field parameters: 



0m (^o) , 

P = ^ (flkta-k;) , 



(10) 

(11) 

(12) 



where (jjm is the molecular field in the closed channel, p is 
the pairing field, and / is half of the number of fermionic 
atoms in the open channel. The interaction Hamiltonian 
'Haa can therefore be written as 



^''3 k^ k 

Similarly, we approximate Ham as 

•Ham ^ -3 XI ('^'""kta-ki 



H.c, 



H.c. 



\p\'~f- 

(13) 



(14) 



Thus, within mean-field the total Hamiltonian is given 
by, 



potential fi. Here, = J2kcr "ktr^kcr + 2 J2k ^k^k is the 
operator of total number of atoms and fi is the chemi- 
cal potential. Using the field operators for Bogoliubov 
quasiparticles [33| . ak+ and ak-, JC takes the diagonal 
form. 



/C = X£;kr4rakr-C/fcs(W' + /')+2(j.-Ai)|0™| 



kr 



E 



(Ck + Ubgf) 



(17) 



where Ck = Ck ~ and the energies of Bogoliubov quasi- 
particles i?k± are given by. 



Cki 



2m 



Ubgf 



(18) 



At temperature T, it is straightforward to write down 
the grand thermodynamic potential. 



" = E 

k 



Ck + Ubgf - 



Ek+ + Ek- 



-UbA\p\' + f 



-2(i.-^)|0™|2-fcBTXln 



kr 



(19) 



The field parameters {(t>m-, P, and /) must satisfy the 
coupled self-consistent equations, dil/df — 0, dVl/dp — 
0, and dVl/d4im = 0. Furthermore, the chemical potential 
is determined by the total number of atoms N, i.e., 



A^^-^^2/ + 20?„ 



(20) 



These four coupled equations can be solved to obtain 
the pairing order parameter A = —{UbgP + .90m) and 
chemical potential pL. 



IV. RESULTS AND DISCUSSIONS 



To clearly contrast the two-channel model with single- 
^'^ ^'^ channel model, we focus on the resonant limit and neglect 

X [([/bgP + 50,„)aj^^aLk; + H.c.j (15) t^^e back-ground interaction. By setting Ubg = (and 
'- therefore A = — (7(/)m)and renormalizing the energy of 

molecules v by using [33|, 

where r = ± is the index of helicity branch. By defining 

an order parameter A = —(UbgP + g(f>m) and rewriting all 2 
the field operators in the helicity basis, the total mean- 21/ — 2vq + — 



field Hamiltonian becomes 



2ek 



(21) 



H = —Ubg (IpP + /^) + 2v\(f}m\'^ + X^ (ekr + Ubgf) h\^^h]^r we obtain the coupled gap equation and number equation 

kr in the two-channel model (r = ±), 

A 



X [e-'^^hi^hl^_^ + e^''-hi_hl^_ + H.c. 



.(16) 



9' 



To determine the variational field parameters (0m, p, and 
/), we diagonalize K, = % — fj,Af by using Bogoliubov 
transformation and calculate the grand thermodynamic 
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where nk± = l/(e^''±/'^'^-^ + 1) is the Fermi-Dirac distri- 
bution function. In contrast, the gap and number equa- 
tions in the single-channel model are given by (l3l |. 



k 



E 



1/2 - n^r 



2E^ 



kr 



1 

2^ 



(24) 



(25) 



respectively. By recalling from Eq. ^ that the effec- 
tive s-wave scattering length in the two-channel model 
is —Anh^as/m = g^/(2zyo), it is clear that the gap and 
number equations in both models have the same struc- 
ture. However, additional terms, — 2/i/5^ and — 2A^/5^, 
appear in the two-channel gap and number equations, re- 
spectively. For a finite SO coupling constant A, if g — > oo, 
l^/g^ and 2A^/g^ go to zero. Then, the equations of 
the two models become exactly the same. This is the 
same as the situation without SO coupling. In other 
words, the single-channel model and two-channel model 
coincide with each other in the broad resonant limit, as 
they should be. However, for a finite channel coupling 
strength g, if A is sufficiently large, deep two-body bound 
state (i.e., rashbon) appears, with a divergent chemical 
potential (/Lt — >■ — c»; see Fig. 1(b) below). Thus, we can 
not neglect 2/x/g^ in Eq. anymore and the two mod- 
els are no longer equivalent. In this strong SO coupling 
limit, we anticipate a qualitative difference between the 
single-channel and two-channel models. 

Let us turn to detailed numerical calculations. For 
simplicity, we consider the case in which the temperature 
is zero and the system is exactly at Feshbach resonance 
(i/Q = and aj^ =0). To characterize the width of Fes- 
hbach resonances, we introduce a dimensionless channel 
coupling constant. 



90 



2m 



1/2 



5, 



(26) 



where kp = (Stt^TV/T/)^/'^ is the Fermi wavelength. We 
take the energy and length in the units of the Fermi en- 
ergy Ep ~ h^kp/2m and the inverse Fermi wavelength 
kp^, respectively. The SO coupling strength is measured 
by the dimensionless parameter X/kp. In experiments, 
the typical magnitude of SO coupling strength is at the 
order of the Fermi wavelength, i.e., A — 0{kp). 

Fig. 1 reports the evolution of the pairing gap (a) 
and chemical potential (b) with decreasing the resonance 
width. For comparison, the prediction of single-channel 
model is also shown by solid lines. As seen from Fig. 1(a), 
for systems with small SO coupling strength, the single- 
channel model and two-channel model give the same re- 
sults when f/o is large enough (i.e., go > 20). In the case 
of ®Li atoms with ~ 600 (the purple dot line in Fig. 
1); we can not see the difference from the single-channel 
prediction. This means that the broad resonance condi- 
tion is always valid for ^Li atoms. However, for a smaller 
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Figure 1: The pairing gap A (a) and chemical potential (b) 
as a function of the Rashba SO coupling strength A for dif- 
ferent atom-molecule coupling go or resonance width at zero 
temperature and at Feshbach resonance. 



(70j e.g., K atoms with ~ 50, the pairing gap deviates 
clearly from the single-channel prediction at the typical 
experimental SO coupling strength \/kp = 3, although 
the two models give essentially the same pairing gap in 
the absence of SO coupling. With increasing the SO cou- 
pling strength, the difference between the two models be- 
comes more significant. For even smaller 170 (i-c, = 5)i 
it is interesting that the dependence of the pairing gap 
on SO coupling strength changes qualitatively. The pair- 
ing gap starts to decrease with increasing SO coupling 
strength and vanishes at sufficiently large SO coupling. 

This dramatic change is somehow not anticipated, as 
the pairing gap is always enhanced by SO coupling in the 
single-channel model. It is closely related to anisotropic 
superfluidity caused by the Rashba SO coupling. As dis- 
cussed in Ref. [Ts'l, due to SO coupling the fermionic 
superfluid has mixed singlet and triplet components. 
The fraction of triplet pairing grows with increasing the 
SO coupling strength. Thus, within the single-channel 
model, the amplitude of pairing gap reflects both sin- 
glet and triplet pairing strengths, and increases as the 
SO coupling increases. In the two-channel model, how- 
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Figure 2: (a) The pairing gap A as a function of the resonance 
width go for different Rashba SO coupling strength X/kp- (b) 
The critical resonance width go (below which the pairing gap 
differs more than 5% from the prediction of the single-channel 
model) for different X/kp- 



Figure 3: The molecular fraction 2\(j}rn\'^/N as a function of 
the Rashba SO coupling X/kp for ''"K atoms where go = 50. 



of the pairing gap A as a function of the resonance width 
go, for some selected values oiX/kp. As go increases, the 
pairing gap A grows rapidly at first, and then saturates 
to the prediction of single-channel model. Quantitatively, 
we may define a critical g^, in such a way that above gc 
the relative difference in the pairing gaps predicted by 
the two models is less than 5%. Fig. 2(b) presents gc as 
a function of SO coupling strength. It gives a qualitative 
phase diagram. Above gc we may safely use the single- 
channel model to describe the Rashba SO coupled atomic 
Fermi gas near Feshbach resonances. While below g^ the 
two-channel model must be adopted. For '*°K atoms with 
go — 50, we find that the single-channel model becomes 
insufficient at a moderate Rashba SO coupling, A ~ 3kp. 



EXPERIMENTAL RELEVANCE 



ever, the most important resonant-interaction Hamilto- 
nian Ham is of s-wave character and hence favors the sin- 
glet pairing. As the triplet pairing is favored by Rashba 
SO coupling, the resonance width and SO coupling have 
opposite effects on the pairing gap and destruct with each 
other. The destruction becomes very pronounced with 
decreasing the resonance width, leading to a completely 
suppressed pairing gap at large SO coupling and narrow 
resonance width. 

The suppression of pairing gap can also be mathemati- 
cally understood from the two-channel gap equation, Eq. 
(P^ . where we may treat 2{vo — /^) as the effective en- 
ergy detuning of Feshbach molecules. By increasing the 
SO coupling, the chemical potential will diverge to — oo, 
and hence the effective energy detuning is pushed up to 
the BCS limit. As a result, the pairing gap becomes sig- 
nificantly suppressed. 

Now, it is clear that the broad resonance limit becomes 
much more difficult to reach in the presence of Rashba 
SO coupling. To quantitatively characterize the broad 
resonance condition, we show in Fig. 2(a) the behavior 



To experimentally test our predictions, we consider 
optical molecular spectroscopy, which projects Feshbach 
molecules onto a vibrational level of an excited molecule. 
The rate of excitations enables a precise measurement 
of the fraction of the closed-channel Feshbach molecules 
in the paired state, although the fraction could be ex- 
tremely small [36|. Near resonance, the paired state may 
be treated as dressed molecules [35, 36], 

Idressed) = e^'^'y^l- Z„, |open) -f yZ^|closed) , (27) 

where Z^. can be identified as the component fraction 
of Feshbach molecules, i.e., Zm = 2|0m|^/A^. As an 
concrete example, in Fig. 3, we show the fraction for 
"'"K atoms (with go ~ 50) as a function of SO coupling 
strength. As the SO coupling increases, the population of 
Feshbach molecule is almost flat at first. However, after 
the coupling reaches a critical value A ~ hp, it grows 
very fast. 

The impact of SO coupling on the population of Fesh- 
bach molecules is best seen in Fig. 4, where we present 
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Figure 4: The molecular fraction as a function of the reso- 
nance width go for different Rashba SO coupling. The inset 
shows Qmax (scc tcxt for definition) as a function of the SO 
coupling strength. 

the fraction as a function of the resonance width at sev- 
eral SO coupling strengths. For a given non-zero SO 
coupling, the fraction is a non-monotonic function of the 
resonance width. By decreasing go from the broad reso- 
nance limit, the fraction first grows then drops to zero, 
as a result of the competition between SO coupling and 
resonance width, as mentioned earlier. In the limit of nar- 
row resonance, the vanishing molecular fraction is con- 
sistent with the suppression in the pairing gap at large 
SO coupling, as shown in Fig. 1(a), due to the relation 
A = —g(f>m,- In contrast, in the absence of SO coupling 
the molecular fraction increases steadily with decreasing 
the resonance width. We may define a characteristic ijmax 
at which the fraction reaches its peak value. As shown in 
the inset of Fig. 4, when X/kp is zero, namely there is no 
SO coupling, f^niax ~ 0, and the population reaches unity 
in the limit of go = 0. As the SO coupling increases, 
5max increases. We emphasize that for small resonance 
width go, even a small SO coupling could lead to a strong 
suppression of the population of Feshbach molecules. 

VI. SUMMARY 

In conclusion, we have investigated a Rashba spin- 
orbit coupled Fermi gas near Feshbach resonances, by 



using a two-channel model. When the spin-orbit cou- 
pling strength is small and Feshbach resonance is broad, 
the two-channel model is equivalent to the single-channel 
model, as we may anticipate [sj, [35|. However, for a 
given resonance width, if the SO coupling strength is 
sufficiently large, these two models are no longer equiv- 
alent. Moreover, for a narrow resonance the pairing gap 
and the fraction of Feshbach molecules are strongly sup- 
pressed by SO coupling. We could test these predictions 
by measuring experimentally the niolecular fraction us- 
ing optical molecular spectroscopy |36J. We have charac- 
terized quantitatively the equivalence of the two models 
by introducing a critical resonance width, above which 
the two models are approximately the same. By cal- 
culating the dependence of the critical resonance width 
on the spin-orbit coupling strength, we have found that 
the single channel model may break down for Rashba 
spin-orbit coupled '^'^K atoms at a moderate spin-orbit 
coupling strength. 

Our results are obtained within mean-field theory, 
which is known to provide a qualitative picture of 
resonantly-interacting atomic Fermi gases. For quanti- 
tative purpose, the crucial pairing fluctuation must be 
included. This may be addressed by using many-body 
T-matrix theories in the future [3a, 123, SOl . 
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